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Fradkin's representation is a general method of attacking problems in quantum field theory, having 
as its basis the functional approach of Schwinger. As a pedagogical illustration of that method, we 
explicitly formulate it for quantum mechanics (field theory in one dimension) and apply it to the 
solution of Schrodinger's equation for the quantum harmonic oscillator. 
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I. INTRODUCTION 

In one of his papersi on functional methods E. S. Fradkin derives a formal solution (Fradkin's representationSi) to a 
class of differential equations which occur in quantum field theory. The solution of the problem reduces to operations 
involving functional derivatives alone. As such, it is complementary to the far more familiar solution in terms of 
Feynman functional integrals. In a different context, Fradkin mentions that for the case of a harmonic oscillator "all 
operations are quite simple and brought to completion." In this pedagogical note we carry out this "simple" task 
(which is perhaps not quite so straightforward as it might appear), as an example of the use of this method, which is 
not widely discussed. Elsewhere, we will apply this method to nontrivial physical examples, particularly to describe 
the physics of magnetic monopoles. 

Fradkin's representation was developed in the context of quantum field theories and in this regard it is an extension 
to Schwinger's functional representation^ for the vacuum to vacuum persistence probability amplitude Z = (0+|0_-). 
The latter object, expressed in terms of external sources, generates all the Green's functions of the theory, and 
therefore defines the theory. Schwinger's expression of field theory in this way was the content of important papers 
published in the 1950si^ 



II. STATEMENT OF THE PROBLEM 



Schrodinger's equation for a quantum mechanical oscillator in one dimension is 



1 



2m dx^ 2 



kx^ 



i>{x,t). 



(1) 



In terms of the parameters uj and a, defined as ^ = a^uj and x ~ TT' '^^^ construct the dimensionless pairs 
of variables. Lot and ax. Without any loss of generality, we can choose uj = \ and a = 1, which renders t and x 
dimensionless. In terms of dimensionless t and x Schrodinger's equation takes the form 
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1 



2dx-^ ^ 



i}{x,t). 



(2) 



The solution to this differential equation for a particular initial condition 

V'(x,0) = 7r-3e-^(^-'^)' 



(3) 



is^ 



l\j{x,t) = 7r-3g-|(^-acost)^g-§[t+2aa:sint-ia2 sin2t]^ 



(4) 



where a is interpreted as the amplitude of the corresponding classical oscillator. The goal here is to reproduce this 
result using Fradkin's representation for Schrodinger's equation. 
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III. FRADKIN'S REPRESENTATION FOR SCHRODINGER'S EQUATION 



In this section we shall keep our treatment general for a time-dependent potential of the form U{x,t). Consider 
the following differential equation, constructed by introducing a suitable source term v{t) for the derivative operator 
in Schrodinger's equation, 



d_ 

dt 



W\xA:v] 



v{t)^-iU{x,t) 



W[x,t-v], 



(5) 



where W[x, t; v\ is a functional of the source v and a fimction of x and t. The solution ^l^{x, t) to the original Schrodinger 
equation is related to W[x, t; v] by the prescription 



ij{x,t)^{W[x,t;v]}^^^. 



(6) 



The initial condition to Eq. ^ is prescribed, in terms of the initial condition to the original Schrodinger equation, to 
be 



W[x,0;v] = ■4'ix,0). 

Notice that we have chosen our initial condition to be independent of v, thus we have 

5 



5v{t) 



W[x,0;v] = 0. 



(7) 



(8) 



Before going ahead, it is good to appreciate the unitary property of the differential equation (0). We know that 
Schrodinger's equation preserves the norm of ip{x,t), 



-[i,,^) = -\dx\^{x,t)\' = Q, 



(9) 



which is the statement of probability conservation. In other words, the time evolution of ip{x^t) is generated by a 
unitary transformation. Using the traditional technique, which involves multiplying Eq. (O by the complex conjugate 
of W[x^ t; v], and then adding the modified equation to its complex conjugated version, we can see that 



dt 



dx 



W[x,t;v] 



0, 



(10) 



if v{t) — v{t)* . (We are also assuming that the potential U is real.) This exposes a very welcome feature of Eq. 
the time evolution involved in W^[a;,t; v] is unitary if v(t) is real. 

To determine Fradkin's representation for Schrodinger's equation we start by integrating Eq. with respect to 
time from to t and obtain the following integral equation 



W[x,t;v] = W[x,0;v]+ f dt' '-^ 



v{t')^^iU{x,t') 



W[x,t';v\ 



(11) 



Using the above integral equation in conjunction with Eq. it is easy to derive a very useful identity 

l"^?^r77 zW[x,t-v]^^W[x,t-v], (12) 

+0 dv(t — e) ox 



which involves the functional derivative of M^[a:,i;f], and e is a positive infinitesimal quantity. This, of course, is the 
reason the source v(t) was introduced. The above identity will play a crucial role in the following development. For 
the sake of dimensional analysis, it might be helpful to remind ourselves that the dimension of a functional derivative 
is not equal to the inverse of the dimension of the function, because the functional derivative is defined as 



5f{x) 
Sfiy) 



= Six - y), 



(13) 



where the delta function has the dimension of the inverse of the variable x. 
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A. Proof of Eq. (IT^ 



We begin the proof of Eq. (|12() by introducing a variable r which is defined in the domain t — e < t < t, where e is 
a positive infinitesimal quantity. Starting from Eq. Ijlll) and using Eq. © we can write 



Sv{t) 



W[x,t;v] 



dt'6(t' -T)-^W[x,t';v] + 
ox 



dt' 



Sv{t) 



W[x,t';v], (14) 



where we have used the definition of functional derivative in Eq. (|13() . The last expression when iterated gives us 
S 



Sv{t) 



W[x,t;v] = / dt'S{t' ~T) — W[x,t';v] 
Jq (jx 



+ 



dt' 



dt"dit" -T)-^W[x,t";v] + 
ox 



(15) 



which is a series expression for the functional derivative of W with respect to v. The first term in the series trivially 
evaluates to our desired result because < {t',T} < t. In evaluating the second term we should be careful because 
the contribution from the delta function depends on whether or not r is greater than t' . We evaluate the second term 
by splitting the integral over t' aX t' = t — e. It is clear that only the integral over t' from t — e to t can contribute, 
and, by the mean value theorem, this integral is of order e. Higher terms in the series are even smaller, so we conclude 
that in the limit e ^ the result (|12|) holds true. 



B. Expression for wavefunction 

To proceed let us write the solution to W[x, t; v] in terms of an auxiliary functional Wi: 

W[x,t;v] = lim 6*^0 Wi[x,t;v], 

e— >0 

which when substituted into Eq. (0 and after the use of Eq. (|12|l gives 



d_ 

dt 



Wi\x,t-v] 



lim e 



v{t)^-iU{x,t) 



(16) 



(17) 



the functional operator having removed the second derivative term in Schrodinger's equation. The exponential of the 
functional operator in the above equation contributes nothing when it acts on the quantity in the brackets due to 
the fact that the sources are defined at different times, that is, because t — e < t. As a consequence, the exponential 
terms in Eq. ((T7|) nullify each other by passing through the bracketed quantity. The differential equation for Wi thus 
takes the form 



d_ 
dt 



Wi[x,t;v] = 



v{t)-^-iUix,t) 
ox 



Wi[x,t; v]. 



Let us eliminate the term by writing 



Wi[x,t;v] = e^odrv{r}£ w^^x^t;v], 



(18) 



(19) 



which when substituted in Eq. H18|l gives 



i^^W2[x,t;v] ^ e-''"'*™(");s[/(2:,t)e-^o'i™(-)* l^2[x,t;i'] 
= U (^x- J^dTv{T),t^W2[x,t;v]. 

The last equation can be immediately exponentiated, 

W2[x,t;v] = W2[x,Q;v] fo dt'u{.^I^' dr vir),t') ^ 



(20) 



(21) 
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Thus, following the sequential substitutions in Eqs. (|16|l . (|19|l . and H21() . we have a formal solution to the differential 
equation in Eq. |SJ|: 



Vl^[a;, t; w] = lim e 



Io<irj^^J^drvir)£^,^ q. dt' u{x~f^' drvir),t') 



(22) 



Using the above formal representation for W[x, t; v] we can write Fradkin's representation for Schrodinger's equation 
for the time-evolved wavefunction using the prescription in Eq. © as [because W{x, 0, v) = Wi{x, 0, v) = W2{x, 0, v) = 
i^{x,0)] 



^Pix, t) = lim "^^^efo '^-(-)i^^(x, 0)6"^ ^« '^*'^(-/o' '^™(-).*') 



t)=0 



(23) 



11=0 



IV. QUANTUM HARMONIC OSCILLATOR 



For the case of the quantum harmonic oscillator we have U{x,t) = 2^ j for which Schrodinger's equation takes 
the form in Eq. (j^J, where we impose the initial condition il^{x,0) given in Eq. 0. Fradkin's representation for the 
quantum harmonic oscillator is obtained by substituting U{x,t) and i/'(a;,0) into Eq. (|23ll . which gives 



For our purpose it will be convenient to rewrite the above equation in the general form 

il}{x, t) = lim TT-^ ^« ^« ^T^IF^-^^*''*") ^wi^e^ ■''0 ■''0 'it"^(t')s(t',t")«(*")+/o' <it',,(t')C(t')+fl(t) 

where the various kernels are 

A{t',t") ^i5{t' ~t") 

B{t',t") = -{l + it<{t',t")} 



D = 



C{t') = -(x- a) - ixt' 



R{t) 



'-{x-aY 



-xH. 



(24) 

(25) 

(26a) 
(26b) 



Here t^{t',t") stands for the function which picks the minimum among t' and t" . In this form, we notice that it 
is possible to carry out the operations of taking the functional derivatives in Eq. H25|) because the exponents are 
quadratic in v. 

To evaluate Eq. (|25ll . consider the discrete analog, constructed from symmetric matrices A and B, column vectors 
c,x and V = and the number r, for which wc have the identity 
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V •A-Vg^x-' -B-x+c^ -x+r 



ix'^BKx+c'^Kx+Ac'^KAc+r+iTrlnK 



where K is the solution to the matrix equation 



1 A B 



K = 1. 



(27) 



(28) 



Generalizing the identity (|27() to functions, and, using it to evaluate the functional derivatives in Eq. H25f) . and at the 
end setting u = 0, we get 



■0(x,t) = TT^J ^R{t) + ^ dt' dt" J* dt"'C(t')K(t' ,t")A{t" ,t"')C(t"') + ^TTlnK ^ 

where K(t',t") is the solution to the integral equation [which is the meaning of Eq. (|28|) for functions] 



dr' 



K{T',t") = d{t' -t"). 



5{t' -t')- / dr" A{t' ,t")B{t" ,t') 
Jo 

Using the expressions for A{t',t") and B{t',t") in Eqs. H26a() and (|26b|l the above equations simplify to 

ij{x,t) ^TT-i e^^*^ + ^ fo dt' lo dt"C{t')Kit',t")C{t") + iTrlnK ^ 

where K{t' , t") and TrlniiT are to be determined using the integral equation 

K{t',t")+i f dT\l+it<{t',T)\K{T,t") =S{t' -t"). 

Jo 



(29) 
(30) 

(31) 
(32) 
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A. Solution of integral equation 



We start by differentiating the above integral equation twice with, respect to t' to obtain the differential equation 
satisfied by K{t',t"): 



dt'' 



K(t'X) 



(33) 

(34a) 
(34b) 

We have presumed that the delta function in Eq. (|32|) does not contribute &i t' — Q and t' = t. (But see below.) In 
terms of a Green's function M{t' , t"), which satisfies 



The boundary conditions on K{t',t") are read off from Eq. (|32|l to be 

K{0,t") = -i f dTK{T,t"), 



K{t,t") = K{S),t")+ / dTTK{T,t") 



we can write 



K{t',t") 



dt' 



dt' 



M{t',t") ^S{t' -t"), 



rM{t' ,t") = S{t' ~ t") ~ M{t' ,t"). 



The continuity conditions for M{t' , t") are dictated by the Green's fmiction equation (|35|) to be 



{M(i',t")},,=,„+, - {M{t'X)}t'^t"-5 - 0, 
{M'(t',t")h,^,„+,-{M'(i',f")},,^,„_, = l, 



(35) 



(36) 



(37a) 
(37b) 



and the boundary conditions on M{t',t") are prescribed by the boundary conditions on K{t',t") in Eqs. (|34a|l and 
(|34b|l . We start by writing the solution to M{t',t") in the form 

A^d' i"\ - / "(^") + co^*'' < < i" < t, 

^^^^^'^ 7j{t")smt' + ^{t")cost', 0<t" <t' < t, ^'^^> 

in terms of four arbitrary constants. Using the continuity conditions H37|l to determine two of the four constants gives 



K{t', t") = 5{t' - t") ~ a{t") sint' - $(<") cost' - sint> cosi<, 



(39) 



where we have suppressed the t' and t" dependence in t^{t' , t") and t>(t', t"). Using the above expression for K{t' , t") 
in Eqs. H34a|l and Ij34b(l we get the equations determining a{t") and £,{t") to be 



which easily yields 



a(t")«[l - cost] + ^(t")[l + « sint] = i cost cos t" - sint", 
Q;(t") cost — ^(t") sint — —cost cost", 



a{t") = -e-"cos(t-t"). 



le 



cost. 



Using Eqs. H41a|l and 141b|l in eqn. H39(l we can obtain the solution to K{t',t") in the form 

K{t', t") = 5{t' - t") - i cos(t - t') cos(t - t") - sin(t - t<) cos(t - t>), 
which can be verified to satisfy the original integral equation H32|) by substitution. 



(40a) 
(40b) 



(41a) 
(41b) 



(42) 



6 



B. Evaluation of Tr In 



We start by observing that K(t',t") and M{t',t") are related by Eq. and M{t',t") satisfies Eq. These 

observations guide us to construct the eigenvalue equation 



dt'' 



+ 1 



(43) 



where we presume that there exists a Hilbert space where ipn{t') are orthonormal eigenfunctions satisfying the or- 
thonormality condition 



and the completeness relation 



(44a) 



(44b) 



We shall defer the derivation of the boundary conditions on ipn{t')- 

It is convenient to evaluate the trace of an operator in terms of its representation as a kernel. Let us denote 



L = T^pY + 1 and introduce the matrix representation 

L = L(i',t") = ^£;„Vn(i')V'n(0*, 

n 

in terms of which our original Eq. H43f) takes the form 

dt"L{t',t")^n{t")=EnMt')- 
In general, we can write the kernel associated with an arbitrary operator F constructed out of L to be 

FiL) = F{t',f) = Y,F{En)Mt')Mt")*- 

n 

As a particular example we can write the kernel introduced in Eq. (|35|1 to be 

= M EE M{t',t") V«(i')V'n(t")*- 

The trace of a kernel in this representation is well defined by 

TrF(L)= / dt'F{t',t') = f dt'y2FiE^)Mt')Mt'r ^y^FiE^). 
Jo Jo „ 



(45) 



(46) 



(47) 



(48) 



(49) 



where in the last step we have used the orthonormality relation (|44a|l . 

Using the above prescription and using Eqs. (|36|l and (|48|l we can write the kernel representation of the operator 
K{t', t") to be 



(50) 



where in the last step we have used the completeness relation (144 b|) . This implies that the eigenvalues of K are 
(1 — 1/En). Therefore, the trace of the logarithm of the kernel K{t',t") is 



TrlnK = E 1^ 



1 



1 

En 



(51) 
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Now we return to the derivation of the boundary conditions on ■!/)„(t'). We multiply ipn{t") by the integral equation 
for the kernel K{t',t") in Eq. H32I) . integrate with respect to t" , and recognize that acting on ip^, K may be replaced 
by its eigenvalue: 

[l-^^Mt')+i dT[l + ii<(i',r)] (^l--^^^„(r)=Vn(i'), (52) 
which further simplifies into 

i!n{t')^~i{l~En) I dr[l + ii<(t',r)lv„(r). (53) 



From the last equation we read off the boundary conditions on ipnit') to be& 

MQ) = -i(l-S„) / drMr), (54a) 
Jo 

Mt) = V'n(O) + (1 - i;„) / dTTMr)- (54b) 

Jo 

Now in terms of = 1 — £'„, the desired trace takes the form 

Tr InK = ^In ^1 - — ^ ^IndetK where detK = J|^l-— ^ . (55) 
The solution to Eq. is 

l/jnit') = AnCOSy/I^it' + BnS'my/I^t', (56) 

where the coefficients An and Bn are in principle determined by the boundary conditions. For the special case of 
Zn — 0, the solution is 

^n{t') ^ Cit' + C2 (57) 

which leads to the trivial solution ipnit') — 0, because the only solution to ci and C2 allowed by the boundary conditions 
are Ci = and C2 = 0. Thus we conclude that z„ is not an eigenvalue. Then we substitute the solutions (|56|l into 
Eqs. H54a|l and (|54b(l to get the equations relating An and i?„: 

An [1 + iy/z^ sin y/z^t] + Bniy/z^ {I - cos y/z^t) = 0, (58a) 
An sin y/z^it - Bn cos y/z^t = 0, (58b) 

which leads to nontrivial solutions to ipnit') only when the eigenvalues satisfy the characteristic equation P{z) = 0, 
where the determinant of the coefficient matrix in Eq. (|58() is 

P{z) = cos y/zt + iy/z sin y/zt. (59) 

As an illustrative example consider a polynomial Q^^\z) of order TV, which can be written in terms of the roots 
zi, Z2, ■ ■ ■ , ZN of the equation Q^^\z) = as 

g(^)(z) = dNiz - zi){z -^Z2)...{z- zn). (60) 

In this form we can deduce 



QW(0) V ^1/ V Z2j"\ ZN 

Presuming that the nonpolynomial P(z) in Eq. H59() has a similar infinite product representation in terms of its roots, 
we can write 

Thus we have established the simple result 

TrlnK = -ii. (63) 
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V. 



RESULTS AND CONCLUSIONS 



Using the solution for the kernel (|42() we can evaluate the first two terms in the exponent in Eq. H31fl to be 




(64) 



Then when we add iTr InK in Eq. ^ we obtain 



■l{j{x, t) — TT 4e 2 




ia^ sin2t] 



(65) 



which is the required solution Q). 

The reader might now rightfully object that we have honed a mighty machine to crack open a peanut. However, 
it seems to us that we have performed a useful pedagogical purpose here in exposing explicitly functional techniques 
which should have much broader applicability in nontrivial contexts. We will present such applications in a subsequent 
publication. 
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® To derive these results from the boundary conditions on the kernels (I34II . we must include the 5- function terms there: in 
effect 
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K{Q, t") K{Q, t") - 5(0 - t") = -M(0, t"), K(t, t") K(t, t") - 5{t - t") 
Then the eigenfunction construction II5UI I correctly implies the boundary condition II54I I. 



M{t,t"). 



